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The metric of a tidally distorted, nonrotating black hole is presented in a light-cone coordinate 
system that penetrates the event horizon and possesses a clear geometrical meaning. The metric 
is expressed as an expansion in powers of r/IZ <C 1, where r is a measure of distance from the 
black hole and 1Z is the local radius of curvature of the external spacetime; this is assumed to be 
much larger than M, the mass of the black hole. The metric is calculated up to a remainder of 
order (r/IZ) 4 , and it depends on a family of tidal gravitational fields which characterize the hole's 
local environment. The coordinate system allows an easy identification of the event horizon, and 
expressions are derived for its surface gravity and the rates at which the tidal interaction transfers 
mass and angular momentum to the black hole. 



Introduction. The numerical integration of the Ein- 
stein field equations for the simulation of a two-body 
system of black holes is currently the subject of intense 
research activity. This work is largely motivated by on- 
going attempts to measure gravitational waves, an effort 
which requires the theoretical input of detailed predic- 
tions for the expected waves. The numerical simulations 
must proceed from initial data that correctly describe a 
situation of astrophysical interest, a spacelike hypersur- 
face that contains two black holes in bound orbital mo- 
tion, together with the gravitational radiation that was 
emitted by the system in the recent past. To find realis- 
tic initial data has proved a challenge |lj , and it has been 
suggested 0, Q that these could be reliably constructed 
using post-Newtonian (and post-Minkowskian) methods. 

While a post-Newtonian metric is adequate to describe 
the mutual gravity of a weakly-bound two-body system, 
it fails to describe the strong gravity of each black hole. 
Building on techniques developed by previous authors 
0, 0| , Alvi remedied this situation [3] by matching the 
post-Newtonian metric to two perturbed versions of the 
Schwarzschild metric (each black hole being tidally dis- 
torted by the other hole), in a buffer region in which both 
approximations are valid. Alvi's metric is ready to use in 
numerical simulations, and in this Letter we improve on a 
part of Alvi's construction by providing a more accurate 
version of the perturbed Schwarzschild metric. 

We calculate the metric of a tidally distorted, nonrotat- 
ing black hole to a level of accuracy that has never been 
attempted before. The hole is immersed in an external 
spacetime whose local radius of curvature is 1Z; this is as- 
sumed to be much larger than M, the mass of the black 
hole. The tidal distortion is created by the Weyl curva- 
ture of the external spacetime, and the Schwarzschild 
metric acquires a correction that can be expanded in 
powers of TlT 1 . We calculate all correction terms of order 
1Z~ 2 and 1Z~ 3 ; while the former were known previously 
0, and used in Alvi's construction [j| , the latter are 
new. The metric is presented in a light-cone coordinate 
system that penetrates the event horizon and possesses 
a clear geometrical meaning. These coordinates, which 



extend over a neighborhood of the black hole that is large 
relative to M but small relative to TZ, allow an easy iden- 
tification of the perturbed event horizon. As an applica- 
tion we calculate the hole's surface gravity and the rates 
at which the tidal interaction transfers mass and angular 
momentum to the black hole. Because of space limita- 
tions we present our results with very few derivations; a 
forthcoming paper will supply the missing details. 

Advanced coordinates based at a world line. We begin 
with the construction of a light-cone coordinate system 
in a neighborhood of a geodesic world line; the only as- 
sumption placed on the spacetime is that its Ricci ten- 
sor vanishes in the neighborhood. This spacetime will 
later play the role of the external spacetime, and the 
world line will later be replaced by the world tube traced 
by the black hole. The quasi-Cartesian version of the 
advanced coordinates is denoted {v,x a ), and the quasi- 
spherical version (v,r, 9 A ). The null coordinate v labels 
past light cones that are centered on the world line, r is 
an areal radius and an afhne parameter on each cone's 
null generators (with r = on the world line), and the 
angles 9 A — (9, <j>) (A — 2, 3) are constant on each gen- 
erator. The coordinates x a (a = 1,2,3) are defined by 
x a := ril a (9 A ) with £l a = (sin 6 cos (f), sin 6 sin (f>, cos 9) . 

Let 7 denote the geodesic world line, which is param- 
eterized by proper time r; its normalized tangent vector 
is u^. Let e£(r) be a triad of orthonormal vectors (all 
orthogonal to it M ) which are parallel transported on the 
world line. Let a{x, x 1 ) be half the squared geodesic dis- 
tance between x and x' , and let a a >(x, x') denote its gra- 
dient with respect to x'. The advanced coordinates are 
defined as follows. Let x' be the point on the world line 
at which r takes the value v. The set of events x such 
that <j(x,x') — and x is in the past of x 1 defines the 
past light cone of x 1 . We assign the advanced-time coor- 
dinate v to all these events. The spatial coordinates of a 
point x on the cone are then x a — —5 ab e^ (v)a a '(x,x'). 
We introduce r := (S a bX a x b ) 1 ^ 2 , which can be shown to 
be an affine parameter on the null generators of the cone 
v = constant. We also introduce fl a := x a /r, which may 
be parameterized by the two angles 9 A . It can further be 
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TABLE I: Irreducible tidal fields. Each field is identified with 
a sans-serif superscript that specifies its multipole content. A 
field labeled with a "q" is a quadrupole field, and one labeled 
with an "o" is an octupole field. The vector and tensor fields 
are all orthogonal to £l a , and all tensor fields are symmetric 
and tracefree. 

£ q = £ cd fl c n d 

si = (s a c - Q a n c )£ cd n d 

£ ab = 2(5 a c - n a n c )(s b d - n b n d )£ cd + (6 ab - n a n b )S" 

Bl = e avq QFB q c n c 

Bi b = e apq n p B q c (s c b - n c n b ) + ( a ^b) 

£° = £ c deQ C Q d Q C 

£°a = (3 a c — Q. a Q c )£cdeQ d Q e 

£° ab — 2(5 a c — Q a Q c )(5 b d — Q, b Q d )£ cde Q e + (5 ab — Q. a Qb)£° 

Bl = le apq n p B q cd ft c n d 

B° ab = le apq n p B q cd (S c b - Q. c Q b )Cl d + (a^b) 



shown that 9 A are constant on the null generators. 

Tensors on 7 can be decomposed in the tetrad (ti M , e^). 
For example, C a0 h := C P e p u x e v h u p and C abc o := 
C 'fj,\ v p e p e b e^u p are frame components of the Weyl tensor 
evaluated on 7. Similarly, C a ow : = C^^^ejJw^ejuV 
and C a0 bo|c : — C '^xv p;uj^a u e b ui: ' e c are frame components 
of the covariantly differentiated Weyl tensor. The 10 
independent components of the Weyl tensor are conve- 
niently related to two symmetric-tracefree (STF) frame 
tensors 4], £ ab := C a0 bo and B a b ■= \e pq C pq bo, where 
£abc is the permutation symbol (frame indices are freely 
manipulated with 5 a b). The components of the differ- 
entiated Weyl tensor can also be expressed in terms of 
STF frame tensors; these include £ a b and B ab as well 

aS £abc • j{C a QbO\c + CcOaO\b + CbQcO\a) and Babe '■ — 
\{ £ a q CpqW\c + Ec q C pqa o\b + S b PC - 'C pgrc 0| o ) . 

These frame tensors allow us to define length scales 
that characterize the tidal gravitational field acting 
at 7. Let TZ be the local radius of curvature, defined 
such that £ a b ~ B ab ~ TZ~ 2 . Let £ be a scale of spatial 
inhomogeneity, defined such that E a b c ~ B abc ~ TZ £ . 
And let T be a time scale associated with changes in the 
tidal field, defined such that £ ab ~ B ab ~ 7£~ 2 T _1 . In 
the sequel we assume, for simplicity, that TZ ~ £ ~ T 
and use TZ as the fundamental length scale associated 
with the tidal field. 

The procedure to obtain the spacetime's metric in the 
coordinates (v, x a ) or (v, r, 9 A ) is detailed in Ref. Q. The 
metric is expressed as an expansion in powers of r/TZ, 
which is assumed to be small; the coordinate domain is 
therefore limited to a neighborhood of 7. The appearance 
of the metric is simplified by forming irreducible combi- 
nations involving the tidal fields £ ab (v), £ ab c{v), B ab {v), 
B a bc(v), the unit frame vector Q a = x a /r, the permuta- 
tion symbol e abc , and the frame metric 8 ab . These com- 
binations are displayed in Table I. 



The quasi-Cartesian form of the metric is 
9w = -l-r 2 £« + ~r 3 £*- ir 3 £° + 0(4), (1) 
9va = n a -^(£ a -Bl) + K 3 (£ a -Bl) 

-lr 3 (£ a -B° a )+0(4:), (2) 
g ab = 5 ab - RA, - ~r 2 (£l b + ^r 3 (£ ab - &J 

-ir 3 (£ ab -B° ab )+0(4), (3) 

where 0(4) stands for "a remainder that scales as TZ~ 4 ." 
The quasi-spherical form of the metric is obtained by 
performing the coordinate transformation x a = ril a {9 ). 
Its nonvanishing components are g vv given by Eq. IjTjI. 
g vr = 1, and 

9vA = -\r\e\-B\)+\r\£\-B\) 

-l r *(£ A -B A )+0(4), (4) 
9AB = r 2 n AB - l -r\£\ B -B\ B ) + ^{£\ B -B\ B ) 

-lr 5 (£ AB -B° AB )+0(4). (5) 

Here n AB := 5 ab n a A Q b B = diag(l,sin 2 0), £\ := 
£\ B := £ q ab n A n b B , and so on, with n A := dn a /89 A . The 
irreducible tidal fields can all be expanded in spherical 
harmonics, and this decomposition is presented in Table 
II. According to Eqs. JT}, @, and ©, ^g = r 2 sin0[l + 
0(4)], which implies that surfaces of constant v and r 
have an area of 47rr 2 [l + 0(4)]. 

Metric of a tidally distorted black hole. The metric of 
an isolated, nonrotating black hole of mass M can be put 
in the form 

glpdx a dx f3 = -fdv 2 + 2dvdr + r 2 n AB d0 A d9 B , (6) 

where / := \-2Mjr and Q. AB d9 A d9 B = d9 2 +sin 2 9 def 2 . 
The coordinates (v,r, 9 A ) have the same meaning as in 
the preceding section: v labels past light cones, r is an 
areal radius and an affine parameter on each cone's null 
generators, and 9 A are constant on each generator. The 
metric of Eq. JHl approaches the Minkowski metric in the 
asymptotic regime r > 2M. 

A tidally distorted black hole is described by a metric 
that is a perturbed version of Eq. ©. In the asymp- 
totic regime r ^> 2M (keeping r <C TZ, which implies 
that M must be much smaller than TZ) this metric must 
approach the metric of Eqs. (JTJ, Q, and JSJ) instead of 
the Minkowski metric. This observation motivates the 
following ansatz for the perturbed metric: 

gw - -f-r 2 e 1 £" + ^r 3 e 2 £"- ir 3 e 3 £° + 0(4), (7) 
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TABLE II: Spherical-harmonic decomposition of the irreducible tidal fields. The first column lists Y lm (9 A ), the real- valued 
spherical-harmonic functions that are used in the decomposition. The abstract index m describes the dependence of these 
functions on <j>\ for example Y 3 ' 2a is proportional to sin 20. The second and third columns define £^(v), £^(v), B q m (v), B^(v), 
the harmonic components of the tidal fields, in terms of their frame components £ a b{v), £ a bc(v), B a b(v), B a b c {v)- The fourth 
column displays the spherical-harmonic decompositions of the quadrupole and octupole tidal fields. The vectorial spherical 
harmonics are Y A m := DaY m and X 1 ™ := —e a b DbY 1 '", where Da is the covariant derivative operator on the unit two-sphere 
(DaQbc = 0) and eab is the Levi-Civita tensor {e$^ = sin#; indices raised with Q AB , the inverse two-sphere metric). The 
tensorial spherical harmonics are Y A "h ■= [DaDb + + l)^As]V im and X% = -\{e%D B + £§ D A )D c Y l ™ . 
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9vr = 1, (8) 

5«A = -|r 3 (e 4 fi-64^) + ir 4 (e 5 q-65^) 

-ir 4 (e 6 ^-6 6 ^)+0(4), (9) 

g AB = r 2 n AB - -r A (e 7 £ AB - b 7 B\ B ) 

+ ^r 5 {e 8 £ AB -b 8 B AB ) 

-l r 5 (e 9 £ AB -b 9 B AB )+0(i). (10) 

The radial functions e n (r), b n (r) must be determined by 
integrating the Einstein field equations; they are inde- 
pendent of 1Z, they are required to be well-behaved at 
r = 2M, and they must approach unity when r 3> 2M. 
The perturbed metric of Eqs. (|7|)- (fTU)) is expressed in a 
gauge that preserves the geometrical meaning of the co- 
ordinates (v,r, 9 A ); this gauge is both physically mean- 
ingful and technically simple. The gauge freedom is not 
exhausted by adopting the ansatz of Eqs. l(7|l- lfTn|) ; re- 
finements can be made when solving the field equations. 

The metric of Eqs. (f^Zfl — (| 1 0|) is used to compute the Ein- 
stein tensor, which is then expanded in powers of TlT 1 . 
There is no contribution at orders TZ° and 7\L _1 , but we 
examine carefully the contributions at orders 1Z~ 2 and 
1Z~ 3 . At order Tlr 2 the Einstein tensor decouples into 
spherical-harmonic components that can be handled sep- 
arately, and it decouples further into even-parity terms 
proportional to £^ and odd-parity terms proportional to 
B* m . The even-parity sector gives rise to a coupled set of 
ordinary differential equations for {ei, 64,67}, while the 
odd-parity sector determines {64,67}. General solutions 



are easily found, and they are constrained by demanding 
that the radial functions be well-behaved at r = 2M and 
approach unity as r — > 00. This leaves a number of inte- 
gration constants undetermined, and those can be freely 
assigned; this represents a refinement of the gauge, and 
we adjust the constants so that a maximum number of 
radial functions vanish at r = 2M. 

This procedure can be completed at order 1Z~ 3 . Here 
the octupole terms (proportional to £^ and B° m ) decouple 
from the quadrupole terms (proportional to £^ and B q m ), 
and as before the even-parity terms decouple from the 
odd-parity terms. We thus obtain coupled sets of equa- 
tions for {e 2 ,e 5 ,e 8 }, {65,^}, {e3,e 6 ,e 9 }, and {6 6 ,6 9 }. 
These are integrated using the same techniques as for 
the TZ~ 2 part of the metric. The radial functions are 
listed in Table III. 

The quasi-Cartesian form of the perturbed metric is 
obtained by inserting e„ and 6„ at appropriate places 
in Eqs. CEJ-©; f° r example, the coefficient of £^ in 
g va becomes — |r 2 e4. The complete determination of 
the metric requires the specification of the tidal fields 
£ab(v), £ a bc(v), B ab (v), and B a bc(v). These characterize 
the black hole's local environment, and they are deter- 
mined by matching the hole's local metric to a metric 
that describes a larger portion of the spacetime 0,0. For 
example, the global metric could be the post-Newtonian 
metric of a two-body system, and Alvi [2| has shown how 
to carry out the matching procedure in this context. 

The perturbed horizon. To locate the event horizon of 
the tidally distorted black hole we note first that the vec- 
tor r a := d a r is normal to any surface of constant r, and 
that it is null when r = 2M : g al3 r a rp = 0(4). The gener- 
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TABLE III: Radial functions that appear in the metrics of 
Eqs. (T71>— (11UI . expressed in terms of x := r/(2M) and / := 
1 - 2M/r. At r = 2M we have e 7 = |, e 9 = ±, b 7 = -§, 
and feg = —jq, with all other functions vanishing. 
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however be calculated using the techniques introduced 
in Ref. ||. Such a calculation reveals that the averaged 
rate at which the black hole increases its mass because 
of tidal heating is given by 

16 / 16 \ 

+ 4725 M8 \ £abc£abC + y SahcSah / (12) 

with a remainder of order (M/1Z) 9 . The rate at which 
the horizon's surface area grows is then (dA/dv) = 
(8tt / k) (dM / dv) . The tidal interaction also transfers an- 
gular momentum to the black hole; the averaged rate of 
change of its angular-momentum vector is given by 

- ^ E M*e a be (e b de S cdB + fB b d eB cde ) (13) 



ators of this null surface have r a as their tangent vector, 
and a simple calculation reveals that r a — t a + 0(4), 
where t a (d/dx a ) := d/dv. The generators are therefore 
parameterized by v and they move with constant values 
of 6 A (up to corrections of order 1Z~ 4 ). That r = 2M 
is in fact the event horizon is confirmed by a theorem 
established in Ref. |8(: Let g a p be the metric of a per- 
turbed Schwarzs child black hole, let h a p '■= g a p ~ g^p 
be the metric perturbation, and let t a be tangent to the 
null generators of the perturbed horizon; in a coordinate 
system in which h a pt@ = on the horizon, the horizon 
has the same coordinate description in the perturbed and 
unperturbed spacetimes. In other words, the gauge condi- 
tions h av = at r = 2M guarantee that the event horizon 
is located at r = 2M in both spacetimes. Because the 
metric of Eqs. lf7|l- (|lU[) satisfies these conditions to order 
lZr A , we conclude that the surface r — 2M[1 + 0(4)] is 
the event horizon of the tidally distorted black hole. 
The surface gravity n of the perturbed horizon is de- 



fined by the statement t a .pt" = nt a , which holds on 
r = 2M and expresses the fact that the integral curves of 
t a are null geodesies. A straightforward calculation gives 



1 

AM 



16 

1 + — M 3 e ab (v)n a n b + 0(4) 



(ii) 



The expansion of the congruence of null generators is cal- 
culated as 8 = {g aP +t a l 13 + l a t' 3 )t a] p, where l a := -d a v 
is the other null normal to the event horizon. A straight- 
forward calculation reveals that 9 = 0(M 4 /1Z 5 ). This 
result implies that within the accuracy of the calculation, 
each cross section v — constant of the event horizon is in 
fact an apparent horizon. 

The slow growth of the event horizon is not revealed by 
a direct examination of the horizon's geometry. It can, 



with a remainder of order M 9 /1Z S . 

Conclusion. The metric of an arbitrary vacuum space- 
time was expressed in a light-cone coordinate system cen- 
tered on a geodesic world line; the metric was expanded 
to third order in powers of r/TZ, which is assumed to be 
small in the coordinate domain. This provided asymp- 
totic conditions on the metric of a tidally distorted, non- 
rotating black hole, which was obtained by integrating 
the equations of black-hole perturbation theory to order 
1Z~ 3 . The coordinates allow an easy identification of the 
event horizon, and the effects of the tidal interaction on 
its physical properties were explored. This work was sup- 
ported by the Natural Sciences and Engineering Research 
Council of Canada. 
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